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Based on suitable video recordings of interactive pedestrian motion and improved track- 
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parameter specifications for the social force model. The calibrated model is then used 
for large-scale pedestrian simulations of evacuation scenarios, pilgrimage, and urban en- 
vironments. 
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1. Introduction 

In recent years, pedestrian research has become a very active field. Queueing- 
theoretical approaches (see, e.g., Refs. [23,25]) and fluid-dynamic models [6,16,17] 
have faciUtatcd insights into the large-scale dynamics of pedestrian crowds, in par- 
ticular into the phenomon of intermittent flows and stop-and-go traffic emerg- 
ing at extreme densities [11]. However, most recent work is dedicated to micro- 
scopic simulation models such as the social force model [8-10,14,38], cellular au- 
tomata [20,21,30], particle hopping models [27], and others, e.g. multi-agent ap- 
proaches [1,37]. 

While some of these models are aiming at the simulation of real evacuation sce- 
narios, others are focussing on an understanding of the self-organized dynamical 
patterns that a realistic model must also be able to reproduce. But how to calibrate 
these models? One common method is to modify model parameters and model com- 
ponents long enough until the simulations fit the fundamental fiow-density diagram 
such as the one proposed by Weidmann [35]: 



V{p) = 1.34ms-^ |l - exp [-1.913m-^ - | 



(1) 
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This method can certainly be appHed, but it does not necessarily lead to realis- 
tic and reliable simulation results. First of all, the fundamental diagram depends 
enormously on the body size distribution [31] and other characteristic features of 
the crowd Second, most models fitting the fundamental diagram must be separately 
calibrated to many differen situations, as they do not reproduce the various dynam- 
ical phenomena in pedestrian crowds well. These phenomena include the formation 
of lanes of uniform walking direction in crowds with oppositely moving pedestrians, 
stripe formation in two intersecting pedestrian streams, and oscillations of the pass- 
ing direction at bottlenecks under normal conditions with moderate densities [9, 12]. 

Apart from the latter two phenomena, which have been discovered only recently, 
the social force model, despite its simplicity, has been shown to reproduce the dy- 
namics in crowds qualitatively well. However, not much effort has been spent in the 
past to calibrate it for the application to dimensioning and evacuation problems. We 
will, therefore, introduce a method to calibrate microscopic pedestrian simulation 
models using trajectory data reflecting pedestrian interactions, and apply it to the 
social force model. 

Our paper is organized as follows: Section [2] will shortly introduce the social force 
model and some specifications of it. In Sec. |3] we will then introduce our method 
to optimize parameters based on trajectory data. For this, we will shortly describe 
the video tracking method applied to gain the trajectory data. Moreover, we will 
discuss how the data can be used to determine the interaction forces themselves. 
Afterwards, in Section [4] we will present some examples of large-scale pedestrian 
and evacuation simulations. Finally, Sec. [5]will summarize and conclude our paper. 

2. Specifications of the Social Force Model 

The social force model can be underpinned with a social science model of behavioral 
changes proposed by Lewin [22] . He has introduced the idea that behavioral changes 
were guided by so-called social fields or social forces. This idea has been put into 
mathematical form by Helbing [4] and applied to opinion formation [7] , pedestrian 
motion [5], and vehicle traffic [33]. The social force model for pedestrians assumes 
that each individual a is trying to move in a desired direction Cq with a desired 
speed and that it adapts the actual velocity Va to the desired one, t^^ — v^e^ 
within a certain relaxation time r^. The velocity Va{t) = dfa/dt, i.e. the temporal 
change of the location Tq, (t) , is itself assumed to change according to the acceleration 
equation 
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The terms faf}{t), fai{t) denote the repulsive forces describing the attempts to keep 
a certain safety distance to other pedestrians (3 and obstacles i. The fluctuation term 
S,a{t) reflects random behavioral variations arising from deliberate or accidental 
deviations from the average way of motion. The above equation are nonlinearly 
coupled Langevin equations and can be solved numerically using Euler's method. 
In very crowded situations, additional physical contact forces come into play [10]. 
For the time being, we will assume a simplified interaction force of the form 

/a/3(i) = w{ipa,[3{t))g{do,/3{t)) , (4) 

where dafj — ra ~ is the distance vector pointing from pedestrian (3 to a and 

ifap the angle between the normalized distance vector dap — dap/\\dap\\ and the 

direction Cq, of motion of pedestrian a, i.e. cos^Lpap) — Bq • d^p. It has, for example, 
been suggested to reflect by the function 

-(^..(t))=(A. + (l-A.)^±i^^) (5) 

that the reaction of pedestrians to what happens in front of them is much stronger 
than to what happens behind them. Here, Aq with < < 1 is a parameter 
which grows with the strength of interactions from behind. Later on, we will try to 
determine this angular dependence from video tracking data. The distance depen- 
dence g{dap{tfj has been specified in different ways. We will start with the circular 
specification of the interaction force, 

g{d^p) = A„e(^°+«^-'^°'')/^°j;;3 , (6) 

where Ra and Rp denote something like the radii of pedestrians a and /?, while 
Ac and are parameters. reflects the strength of interaction, while Ba cor- 
responds to the interaction range. While the dependence on a explicitly allows for 
a dependence of these parameters on the single individual, we will assume Aa = A 
and Ba = B va the following. Otherwise, it would be hard to collect enough data 
for parameter calibration. 

3. Evolutionary Model Specification Via Video Tracking 

Let us now turn to the question how the parameters A and B and the angular 
interaction strength w^ipap) can be determined from video recordings. 

3.1. Video Tracking Method 

We have made several video recordings of pedestrian crowds in different natural 
environments in Budapest (Hungary) and Stuttgart (Germany) from the top (see 
Fig. [ij. The dimensions of the recorded areas was known, and the floor tiling or 
environment provided something like a "coordinate system". The heads were au- 
tomatically determined by seaching for round moving structures, and the accuracy 
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Fig. 1. Video recordings used for parameter calibration in this study. Top: Entrance and exit area 
of two escalators. Middle; Free movement in a shopping mall in Budapest. Bottom: Pedestrian 
crossing experiment with students (see Ref. [9] for details). At least in the first two cases, pedes- 
trians did not know they were recorded. Therefore, it can be assumed that their behavior was not 
influenced and representative for middle European conditions. 



of tracking was improved by comparing actual with linearly extrapolated positions 
(so it would not happen so easily that the algorithm interchanged or "lost" closeby 
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pedestrians). The trajectories of the heads were then projected on two-dimensional 
space in a way correcting for distortion by the camera perspective. A representative 
plot of the resulting trajectories is shown in Fig. [2j However, it should be noted 
that extracting trajectory data from pedestrians is nothing new. It has been done 
in the past with infra-red sensors [19] or video recordings [15,31]. 




Fig. 2. Video tracking used to extract the trajectories of pedestrians from video recordings close to 
two escalators. Left: Illustration of the tracking of pedestrian heads. Right: Resulting trajectories 
after being transformed onto the two-dimensional plane. 



3.2. Parameter Optimization by an Evolutionary Algorithm 

For model calibration, we have used a hybrid method fusing empirical trajectory 
data and microscopic simulation data of pedestrian movement in space. To each 
tracked pedestrian, we have assigned a virtual pedestrian in the simulation domain. 
We have then started a simulation for T = 1.5 seconds, in which one pedestrian 
a was moved according to a simulation of the social force model, while the others 
were moved exactly according to the trajectories extracted from the videos. This 
procedure was performed for all pedestrians a and for several different starting 
times t, using a fixed parameter set for the social force model. 

Each simulation run was performed according to the following scheme: 

(1) Define a starting point and calculate the state (position r^, velocity Va, and 
acceleration — dva/dt) for each pedestrian a. 

(2) Assign a desired speed w° to each pedestrian. In our simulations, we have spec- 
ified it by the maximum speed Va{t) = ||-!TQ,(t)|| during the pedestrian's tracking 
time, which is sufficiently accurate if the overall pedestrian density is not too 
high and the desired speed is constant in time. 

(3) Assign a desired goal point for each pedestrian. We have assumed is would 
correspond to the point at the end of the trajectory. 

(4) Given the tracked motion of the surrounding pedestrians /3, simulate the trajec- 
tory of pedestrian a based on the social force model during T seconds of time. 
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starting at the actual location ra{t). 
After each run, we determined the relative distance error 



(7) 



ll^racked^^ _j_ ^racked^^j|| 

After averaging the relative distance errors over the pedestrians a and starting times 
t, the negative value of the result was taken as fitness of the parameter set used 
in the pedestrian simulation^ Hence, the best possible value of the fitness was 0, 
but any deviation from the real pedestrian trajectories would imply negative fitness 
values. 

Moreover, with the fitness measure outlined above, we have used a simple evolu- 
tionary algorithm to obtain the parameter set with the highest fitness value, which 
was —0.39. This value reflects in some sense the stochasticity and/or heterogeneity 
of individual behaviors. For comparison, we also made simulations with A — 0, i.e. 
without any social forces, assuming that persons will never change their velocity. 
This measurement resulted in a fitness value of —0.66. 



3.3. Improved Specifications of the Social Force Model 
3.3.1. Velocity- Dependent Interaction Forces 

Interactions among pedestrians are actually more complicated than suggested 
above. For example, it is know that the angle Sap matters, at which two pedes- 
trians a and (3 approach each other (which is given by cos(5q/3) = • e*/?)- Apart 
from this, the step sizes and, therefore, the speeds matter as well. In the following, 
we will shortly discuss two anisotropic models of pedestrian interactions: 

Elliptical specification I: In Ref. [13], a generalization of Eq. ^ was formu- 
lated, which assumed that the repulsive potential 

Vo.p{b) = ABc-^-f''^ (8) 

is an exponentially decreasing function of b with equipotential lines having the form 
of an ellipse directed into the direction of motion as shown in Fig. [3] The semi-minor 
axis hap was determined by 

2&a/3 = \j{\\dai3\\ + \\dap-vpMep\\Y -{vfi At)^ (9) 

in order to take into account the length vp At of the stride (step size) of pedestrian 
P, where vp = \\vf3\\. The reason for this specification was that pedestrians require 
space for movement, which is taken into account by other pedestrians. 
The repulsive force is related to the repulsive potential via 

gc.p{d:ap) = -y,-Vap{bap) = -^^^^^v,- 6„^((i;^) . (10) 

''Note that we have actuaUy averaged over the 30% central values of the relative distance errors only 
in order to remove outliers due to untypical, unreasonable, or exceptional pedestrian behaviors. 
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Considering the chain rule, ||z| 
expHcit formula 



2^, and Vi-||zj| = z/\/z^ — z, this leads to the 



with ijap = vp Ate p. For At = 0, we regain the expression of Eq. (|6|. 



(11) 



b 




Fig. 3. Illustration of the elliptical specification of pedestrian interaction forces. 



Elliptical specification II: Recently, a variant of this approach has been pro- 
posed [32], assuming 



2b ^Ji\Ki,\\ + \Kp - (vp - z7„)At||)2 - [{v/, - v^)At]^ . (12) 

The special feature of this approach is its symmetrical treatment of both pedestrians 
a and p. 

In order to compare the elliptical specifications with the circular one, let us make 
the following thought experiments: 

• Scenario 1: Consider a pedestrian walking towards a wall with a given speed 
but four different orientations (with k £ {1,2,3,4}), as shown in Fig.|4j 
This scenario can be treated analogously to the movement relative to a standing 
pedestrian /?, which implies vp = 0. Then, in case II the repulsive force is a 
monotonously decreasing function of ba/3 given by 



2&a/3 = \Ji\Kt}\\ + ||da/3 + i;„Aie„||)2 - (w„At)2 (13) 



according to Eq. ( 12 ). For all four cases, the values of d :~ ||da/3|! and VaAt (the 
step size of pedestrian a) are the same, but the values of d'^ :~ Wdap + VaAte'^W 
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Fig. 4. Illustration of scenario 1, in which pedestrians are walking at a given speed in different 
directions close to a wall. 



are difJerent. We have di < dl ^ < dl, so that we find > Fj^ = > F^ 
for the magnitudes of the repulsive forces triggered by the wall (as the direction 
of the forces is perpendicular to the wall in all four cases.) This agrees well 
with experience, i.e. the anisotropic and orientation behavior of pedestrians are 
realistically reproduced by the elliptical force specification II. In contrast, the 
elliptical model I implies b^fj = d^p according to Eq. (j9]) and predicts the same 
force in all four cases, as does the circular model. 

Scenario 2: Consider single pedestrians a walking towards a wall, but with 
2 different velocities (u^ < u^) as shown in Fig. [sj Obviously, the elliptical 
specification II predicts that the pedestrian a with the higher speed would 
decelerate earlier, as expected. In contrast, the elliptical model I and the circular 
interaction model are insensitive to the own speed Va- 



T 



«B 



\ 



Case 1 



Fig. 5. Illustration of scenario 2, in which pedestrians are walking with different speeds towards a 
wall. 
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• Scenario 3: Consider the case of two pedestrians walking towards each other 
with different velocities Va > vp, as shown in Fig.|6] While the circular interac- 

a R 

Fig. 6. Illustration of scenario 3, in which pedestrian are walking towards each other with different 
velocities. 



tion force does not imply any velocity-dependence of pedestrian behavior, the 
elliptical formulations do, which is certainly more realistic. The first elliptical 
model assumes a velocity-dependence of pedestrian interactions, but no sensi- 
tivity with respect to the own speed. Therefore, specification II is expected to 
be more realistic. 

For the above reasons, we will not consider the elliptical specification I anymore, 
but focus on the elliptical model II. 

Figure [7] represents the resulting fitness values as a function of the choice of the 
interaction strength A and the interaction range B, assuming w^ip) = 1. 

1 

0.8 
0.6 

< 

0.4 
0.2 



"0 1 2 3 4 5 

B 

Fig. 7. Goodness of fit for different parameter values A and B, when pedestrian trajectories are 
simulated with the social force model and the the force specification (j6]l. The darker the shade, 
the better the fitness is, for the corresponding combination of A and B values. 
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Table 1. Optimal parameter values resulting from our evolutionary parameter optimization for 
three different specifications of the interaction forces between pedestrians (see main text) The 
parameter calibration was based on the 3 video recordings mentioned in Fig.^and the fundamental 
diagram of Weidmann, see Eq. llj. 



Model 


A 


B 


A 


Fitness 


Extrapolation 









-0.66 


Circular 


0.42 ± 0.26 


1.65 ± 1.01 


0.12 ± 0.07 


-0.60 


Elliptical I 


0.11 ± 0.01 


1.19 ± 0.45 


0.16 ± 0.04 


-0.59 


Elliptical II 


0.04 ± 0.01 


3.22 ± 0.67 


0.06 ± 0.04 


-0.39 


Circular 


0.11 ± 0.06 


0.84 ± 0.63 


1 


-0.65 


Elliptical I 


1.52 ± 1.65 


0.21 ± 0.08 


1 


-0.67 


Elliptical II 


4.30 ± 3.91 


1.07 ± 1.35 


1 


-0.47 



One result of our parameter optimization was that, for each video, there was 
actually a broad range of parameter combinations of A and B which performed 
almost equally well, while the optimal value of the anisotropy parameter was A = 
0.11 ± 0.07. This allowed us to apply additional goal functions in our optimization. 
We used this fact to determine among the best performing parameter values such 
parameter combinations, which performed well for all three video recordings, using 
a fitness function which weighted the fitness reach in each single video equally (i.e. 
with a factor 1/3). 

The resulting parameters are listed in Table [l] The parameters corresponding to 
the best fitness value, i.e. the elliptical specification II with A — 0.04, B = 3.22, A = 
0.06, were used in a microscopic simulation with periodic foundary conditions, and 
the resulting fundamental diagram, is shown in Fig. [8] The solid line is the Wei- 
dmann curve, Eq. (jlj. Note that, with additional fine tuning, it will be possible 
to make an even better fit to the Weidmann fundamental diagram, but here we 
rather just show that extracting parameters from videos and perform a microscopic 
simulation will correspond to dynamics close to the empirical fiow-density relation. 

We would, however, like to mention that further velocity-dependent specifica- 
tions of pedestrian interaction forces were proposed in the past, for example, by 
Molnar Ref. [8,26], Okazaki [28], Hoogendoorn [14], and Yu [38]. 

Although, from a theoretical perspective, velocity-dependent force specifications 
are more promising than the circular interaction force introduced in the beginning, 
it turns out that the prediction of pedestrian motion does not become a lot better 
by including additional, speed-dependent parameters. This is illustrated by the 
fitness values in Table [T] Obviously, the deviation of the fitness value from is 
mainly caused by the heterogeneity of pedestrian behaviors. A similar conclusion has 
been drawn from studies specifying behavioral parameters in different car-following 
models [3]. 
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Fig. 8. Comparison of the fundamental diagram for pedestrian flows by Weidmann, see Eq. |[TJ 
with the fundamental diagrams resulting in simulations with the social force model, Elliptical II, 
for the parameter combinations listed in Tab. ^ 



3.3.2. Angular Dependence of Pedestrian Interactions 

Our evolutionary fitting method can be also used to determine interaction laws with- 
out prespecified interaction functions. For example, one can determine the angular 
dependence of pedestrian interaction directly. In order to do this, we have chosen an 
approach using a polygon with n edges in the directions ip € 27rj/n, i S [0, n— 1]. The 
distance of the curve in direction ipi from the origin represents the (relative) angular 
interaction strength w{ipi) with a maximum value of 1. To get some smoothness, 
linear interpolation has been applied between each of the edges. Figure [TO] shows 
the angular dependence of the interaction strength determined from pedestrian 
trajectories. Apparently, pedestrians are only sensitive to what happens inside an 
180 degree angle in front of him or her, which roughly corresponds to the visually 
perceived area. 



3.3.3. Distance Dependence of Pedestrian Interactions 

In a similar way, we have determined the distance dependence of pedestrian inter- 
actions without a pre-specified function. For this, we have adjusted the values of 
the force at given distances dk = kdi (with k G {1,2,3,...}) in an evolutionary 
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Fig. 9. Another good feature of the eUiptical model II is that it produces a non-regular density 
distribution, which is comparable to how persons distrubute in dense crowds, in reality. Left: A 
snapshot from a simulation with the elliptical model II. Right: A snapshot from reality, where 
persons are identified from a picture of a dense crowd. In both cases, the persons are walking from 
left to right. 
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Fig. 10. Angular dependence of the influence of other pedestrians. The direction along the positive 
X axis corresponds to the walking direction of the pedestrian, y to the perpendicular direction. 



way. The resulting fit curve is presented in Fig. [TT] According to our results, the 
empirical dependence of the force with distance can be well fitted by an exponential 
decay, as assumed before. Note that similar studies have been carried out for vehicle 
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Fig. 11. Empirically determined distance dependence of the interaction force between pedestrians. 
An exponential decay fits the empirical data quite well. The dashed fit curve is made with Eq. 
|6]and A = 0.53, B = 1.0. From the simulation with the circular model the resulting parameters 
A = 0.42, B = 1.62 are quite close. 

interactions, recently [34]. However, the situation for pedestrians is more difficult 
due to their two-dimensional motion. Therefore, our evolutionary fitting procedure 
seems to be superior to the Fokker-Planck approach (or Random Matrix Theory 
approach), which was used to calibrate vehicle interactions. 

4. Large- Scale Pedestrian and Evacuation Simulations 

In the following section, we will present some examples for large-scale pedestrian 
simulations with the social force model. For simplicity and computational speed, we 
have assumed circular pedestrian interactions, here. Compared to previous imple- 
mentations of the model, numerical efficiency has been increased a lot, for example 
by cutting off interation forces at a certain, large enough distance (around 5 m) 
and book-keeping of neighboring pedestrians. In this way, the simulation of 30,000 
pedestrians or more in real-time is no problem. Larger crowds have been simulated 
off-line or using parallel computing. 

For example, a simulation have been performed for an evacuation of two decks 
of a ship connected by a staircase. The underlying study has shown that the model 
is also suitable for treating different personalities (the different characteristic be- 
haviors in emergency situations), complex geometries, and three-dimensional inter- 
actions [36]. A simulation [14] have been performed for multi-destination flows at 
Schiphol airport in the Netherlands, while [24] have simulated the evacuation of 
an exhibition area, both with the social force model. For a similar simulation see 
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Fig. 12. Snapshot of an evacuation simulation of a convention or exposition hall. 



Our model has also been used to simulate pedestrian flows on the Jamarat 
Bridge, which has to cope with 3 million pilgrims within a single day (AlGadhi and 
Mahmassani, 1990). Figure 13 illustrates the scenario for the old Jamarat bridge. 
It shows the development of zones of high density and pressure behind the previous 
circular basins. In the middle of each of the three basins is a pillar symbolizing 
the devil. These pillars are to be stoned with seven pebbles each, which takes time 
and delays pedestrian motion. This, in turn, can cause extreme crowding. Note, 
however, that the situation around the pillars has significantly improved since the 
circular stoning areas were replaced by elliptical ones. 

As indicated before, the implementation of our pedestrian microsimulator has 
been optimized for efficiency and scalability in terms of computational speed and 
memory. It is no problem to simulate more than 30,000 pedestrians in real-time on 
a normal PC. Parallel computing even allows one to simulate more than 100,000 
pedestrians in real-time. This is suitable for the simulation of large scenarios such 
as pedestrian flows in extended urban areas (see Fig. 14 1 or during mass events such 
as street parades or carnivals [1]. 

In order to be suitable for scenarios of this complexity, we have developed a 
powerful and flexible pedestrian simulator called "UNIVERSE", in which we have 
implemented certain additional features complementing the social force model: 
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Fig. 13. Top: Simulation of the Jamarat bridge which has to cope each year with miUions of pilgrims 
within a few days. The pilgrims need to get close to the circular area in order to throw stones 
at the Jamarah (the pillar in the center representing the devil). This leads to extreme pedestrian 
densities behind the pillars (dark red and light yellow areas) which have led to fatalities in the 
past. Note that the recent replacement of the circular Jamarahs by elliptical ones has lead to 
considerable improvements in the capacity and safety of the stoning process. Bottom: Photograph 
of the (dark) compression area behind a circular Jamarah. 

(1) a powerful scenario generation module based on XML, 

(2) path finding and route choice modules to describe the selection of destinations, 
including spontaneous stops and impluse shopping, 

(3) calibration modules to adapt simulations to real measurements, 

(4) a statistics module evaluating various performance measures (such as densities, 
velocities, flow efficiency, pressure, overall evacuation time, etc.) to assess and 
compare different scenarios, 

(5) an optimization module based on evolutionary algorithms which allows to im- 
prove the geometrical boundary conditions of pedestrian facilities, 

(6) various visualization modes for the illustration of the simulation results, both 
on- and off-line, and 
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Fig. 14. Densities of pedestrian streams in a simulation of the city center of Dresden, Germany. 
Our simulation handles 30.000 pedestrians. 

(7) a functionality reflecting density-dependent herding effects. 
5. Summary and Outlook 

Tracking algorithms have become a powerful tool for empirical data evaluation, 
recently [15,19,31]. In our study, we have further developed and applied this tool 
to determine trajectories of many interacting pedestrians. These trajectories were 
used to calibrate a microscopic pedestrian model, namely the social force model. For 
this, we have simulated single pedestrians, given the measured motion of the other 
pedestrians. The difference between simulated and measured pedestrian behavior 
has then been minimized through evolutionary variation of the model parameters. 
However, it was also possible to determine the distance- and angle-dependencies of 
the interactions directly. In this way, the previously assumed exponentially decay 
of the interaction strength with distance was empirically confirmed. 

It turned out that the simplest specification of the interaction forces in the social 
force model almost performed as well as velocity-dependent specifications with more 
parameters. We have, therefore, used the circular specification for our large-scale 
simulations, which is in favor of computational speed. With some simplifications 
for the sake of numerical performance, it is no problem to simulate scenarios with 
30,000 pedestrians on-line on a normal PC. Larger crowds can be simulated off-line 
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or using parallel computing. This allows one to study pedestrian zones, convention 
halls, shopping malls, airports, railway stations, city festivals, etc. There are, by the 
way, also special algorithms for the evolutionary optimization of highly frequented 
areas [18]. Hence, pedestrian simulations are now a useful tool for the planning of 
pedestrian facilities and mass events, if properly calibrated. 

Note, however, that that measurement and calibration results should be trans- 
fered from one place to another only with care, as the maximum flow and density 
values are largely influenced by the body size distribution. Other factors may play 
a role as well, for example the age distribution, the fraction of handicapped people, 
the carrying of luggage, alcohol consumption, cultural habits, use of mobile phones, 
and the situational context [29] (e.g. business, leisure, sports, or religious activities). 
All of this should be taken into account in the dimensioning of pedestrian facilities, 
safety (evacuation) analyses, and the calibration of computer models. 
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